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ABSTRACT 

We report the clearest detection to date of dusty torus signatures in a Weak-Line 
Radio Galaxy (WLRG). The deep Spitzer InfraRed Spectrograph (IRS) rest-frame 
mid-infrared (MIR) spectrum of the WLRG PKS 0043-42 (z=0.116) shows a clear 
spectral turnover at A > 20 pm suggestive of warm dust, as well as a 9.7 urn sili- 
cate absorption feature. In addition, the hard X-ray results, based on Chandra data, 
strongly support a picture in which PKS 0043—42 has a torus and accretion disc 
more typical of Strong-Line Radio Galaxies (SLRGs). The MIR and X-ray spectra 
are markedly different from those of other WLRGs at similar redshifts, and here we 
show that the former can be successfully fitted with clumpy torus models with param- 
eters characteristic of Type-2 AGN tori: close to edge-on (i=74°) and relatively broad 
(o-=60°), with an outer radius of 2 pc, N#=1.6±[5; 2 xl0 23 cm" 2 , and AGN bolomet- 
ric luminosity = 1.6±o' 2 xl0 44 erg s" 1 . The presence of a compact torus in 

PKS 0043-42 provides evidence that this WLRG is fuelled by cold, rather than hot, 
gas accretion. We suggest that WLRGs are a diverse population, and PKS 0043-42 
may represent a type of radio galaxy in which the AGN activity has been recently re- 
triggered as a consequence of intermittent gas supply, or in which the covering factor 
of the Narrow-Line Region (NLR) clouds is relatively low. 

Key words: galaxies: active - galaxies: nuclei - galaxies: spectroscopy - galaxies: 
individual (PKS 0043-42). 



1 INTRODUCTION 

Most powerful radio galaxies (PRCs) can be classified 
as Narrow- and Broad-Line Radio Galaxies/quasars (NL- 
RGs, BLRGs/QSOs) on the basis of their spectral fea- 
tures. The unified model for active gal actic nuclei (AGN; 
lAntonuccilll993l ; lUrrv fe Padovanilll995l b proposes the ex- 
istence of a pc-scale obscuring toroidal structure to ac- 
count for the observed differences between their spectra. 
However, there is also a third class, the WLRGs (also 
known as Low- Excitation Galaxies; LEGs), whose optical 
spectra are dominated by the stellar continua of the host 
galaxies (|Laing et al.lll994l : iTadhunter et al.lll998l ') without 
the prominent emission lines characteristic of NLRGs and 
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BLRGs (which can be grouped as SLRGs, or alternatively, 
High- Excitation Galaxies; HEGs). 

The nature of the physical parameters which lead to the 
division between SLRGs and WLRGs still remains unclear. 
Previous studies failed to explain the differences between 
the t wo types as due to differen t or time-varying accretion 
rates (Ghiscllini fc Celotdl200lh . An alternative interpreta- 
tion is that SLRGs are powered by cold gas accretion, while 
WLRGs are fuelled by accretion of hot g as provided by the 
reservoir of their X-ray gaseous coronae (I Allen et al.ll2006l: 
Best et alJ l200rj; lHardcastle et all 120071 ; iBalmaverde et all 
20081 : iButtiglione et alJl201Ch . The high temperatures of the 
hot gas would prevent the formation of the "cold" structures 
(e.g. the Broad-Line Region and the torus). This hypothe- 
sis explains the fact that WLRGs do not normally show the 
high l evel of X-ray absorptio n expected from a molecular 
torus jHardcastle et al.ll2009T ). 
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SLRGs are almost invariably associated with Fanaroff- 
Riley II sources (FRII). On the other hand, while most 
WLRGs show Fanaroff- Riley I (FRI) radio mo rphologies, 
some are classified as FRIIs (|Laing et al.| [l994). Evidence 
for nuclear warm du st emission has been found in the 
MIR for FRII/SLRGs l|Ogle et alJIiooi Ivan der Wolk et all 
120101 ) and, to a lesser extent, for a few FRI and FR II WL- 
RGs (|Leipski et all 120091 ; Ivan der Wolk et al.l l201oh . How- 
ever, none of the WLRGs so far studied have provided 
definitive evidence for a warm compact obscuring region 
e.g. in the form of a spectral turnover at ~20 u.m and a 
clear silicate absorption feature. This lack of conclusive de- 
tection of torus signatures in previous studies may be due 
to contamination of the MIR data by stellar photospheric 
and starburst-heated dust components, which generally pro- 
duces a broad minimum in the 8-10 \im range and can be 
confused with a silicate absorption feature. Strong PAH (in 
particular the 8.2 and 11.3 [im features) and non-thermal 
emission from the synchrotron-emitting core sources can 
also mask the torus signatures. Based on extrapolation of 
the radio core data, the latter emission appear s particularly 
important in WLRGs at MIR wavelengths dDicken et al.l 
l2008l : lvan der Wolk et al.ll2010l ). All of these diluting compo- 
nents need to be accurately subtracted before the presence 
of warm dust can be deduced. 



At optical wavelengths the galaxy PK S 0043-42 is clas- 
sified as a WLRG at redshift z=0.116 l)Tadhunter et all 
Il993h . whereas at radio w avelengths it s hows a clear 
FRII morphology (iMorganti et ai1ll993l , ll999l ). According to 
iTadhunter et all l| 19981 ). the WLRGs are defined as a group 
of galaxies with [O III]A5007 emission line equivalent widths 
below IOA. PKS 0043-42 not only fullfills this criterion, but 
also has an [O III]A5007 luminosity a factor of 10 lower than 
those of the general population of SLRGs a t similar redshifts 
and radio powers (Lfo/// i=5xl0 40 erg s" 1 ; ITadhunter et al.l 
1 19981 ; [Dickon et al] [2009) , as well as line ratios indicating 
a low ionization state (e.g., [O III1A5007 /[O IL1A3727 < 1; 
ITadhunter et al.lll998l ; iLewis et al]|2003l ). The galaxy is at 
the center of a gr oup or cluster, and from NIR imaging, 
llnskip et all (|2010h reported a central isophotal twist and an 
excess of emission along the NS direction. They claimed that 
both are likely to be associated either with the presence of a 
dust lane, or with the apparent interaction with the compan- 
ion object ~15 kp c to the North. However, the deep optical 
image presented in lRamois Almeida et~aI1 l|2010h does not re- 
veal any dust features in PKS 0043-42, but appears instead 
to confi rm the presence of the brid ge with the companion 
galaxy (|Ramos Almeida et"ai]|2010h . 



In this study we report the detection of torus signatures 
in the MIR and X-ray spectra of the WLRG PKS 0043- 
42. In the following sections we describe the observations, 
the SED modelling, the main results and their implications. 
Throughhout this paper we assume a cosmology with Ho=73 
km s _1 Mpc -1 ; f2 m =0.27, and Qa~0-7S, which provides a 
luminosity distance of 516 Mpc for PKS 0043-42 (from the 
NASA/IPAC Extragalactic Database; NED). 



2 OBSERVATIONS 
2.1 Spitzer observations 

MIR spectroscopic observations of PKS 0043-42 were ob- 
tained using the IRS instrument (|Houck et al.l l2004h on 
the Spitzer satellite on 2008 Dec 12 under the program 
2JYSPEC/50558 (PI: C. Tadhunter). The observations were 
made in staring mode using the two low-resolution (R~60- 
120) IRS modules: the short-low (SL; 5.2-14.5 ^.m) and 
the long-low (LL; 14-38 urn). The slit widths are 3.7" and 
10.7" for the SL and LL modules respectively. Four cy- 
cles of 60 s exposures were taken in each of the two SL 
modules, and two cycles of 120 s for the LL modules. The 
data were downloaded from the Spitzer archive in basic cali- 
brated data (BCD) product form, processed with the S18.7.0 
pipeline, and reduced and extracted using the SMART v. 8. 1.2 
program developed by the IRS Team at Cor nell University 
l|Higdon et al] 12004 iLebouteiller et al.ll2010l ). In extracting 
the fluxes we used the optimal extraction function which 
uses a super sampled PSF and weights t he extracted spec- 
tra u sing the signal-to-noise of each pixel (|Lebouteiller et all 
120101) . Note that the IRS pipeline automatically accounts for 
variable slit losses with wavelength. For each cycle the nod 
positions were subtracted in order to perform background 
subtraction on the image, as well as to counteract intrinsic 
sensitivity variations across the detector. For a more de- 
tailed description of the observations and data reduction 
see Dicken et al., in prep. The two resulting spectra were 
then combined to produce the final spectrum shown in Fig- 
ure [T] The absolute uncertainty in the flux calibration at 
the la level is ~15% for the SL and LL data, whereas for 
the spectral shape it is much less than that (|Decin et al.l 
|2004 iRoellig et al]|2004l '). We have scaled the observed IRS 
spectrum to the Multiband Imagin g Photometer for Sp itzer 
(MIPS) 24 urn flux (11.1±0.2 mJv: lDicken et al.ll2008T l. ob- 
tained using a 15" aperture. The discrepancy between the 
IRS and MIPS fluxes at 24 y.m is only ~10%. 

In addition to the IRS spectrum, we compiled: 

(1) A NIR upper limit flux of 0.45±0.02 mjy derived from 
the GALFIT analy sis of Ks-band Sofl images presented in 
llnskip et al. (2010). This value corresponds to an unresolved 
component of 0.93" of FWHM. 

(2) A MIPS 70 urn flux of 9.9±3.0 mjy from lDicken et all 
(2008). For the analysis performed here, we consider this 
flux an upper limit, because at this wavelen gth and with 
the a perture employed for calculating it (24": iDicken et all 
2008) the emission is likely to be dominated by dust heated 
by stars and/or AGN heated dust in the NLR. In Figure 
[2] we present the MIPS images of PKS 0043-42 at 24 and 
70 (im (from the 2Jy web pagc: http://2jy.extragalactic.info 
and described in lDicken et al . 2008), to show the diference 
in the signal-to-noise level of the two detections. 

(3) High angular resolution radio core fluxe s of 4. 5=1=0.8 
mjy at 2.3 GHz (13 cm) from iMorgant i et all (Il997l). an d 
1.4±0.3 mjy at ~15 GHz (2 cm) from lDicken et all (2008). 



The IR-to-radio SED of PKS 0043-42 is shown in Fig- 
ure|3] From the radio core data points it becomes clear that, 
for any reasonable extrapolation of the radio core spectrum, 
non-thermal emission does not contribute significantly to 
the MIR emission of this gal axy. By looking at the h igh res- 
olution radio maps shown in lMorganti et all (|l997l ) we can 
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Figure 1. Low-resolution IRS rest-frame spectrum of PKS 0043- 
42 in relative F„ units (top). Note the emission bump peaking 
at ~20 urn and the 9.7 urn silicate feature in absorption (short 
vertical line), both characteristic of a dusty torus. The spectral 
features and shape are different to those of the 2 Jy median WLRG 
spectrum (bottom). Typical AGN emission lines are labelled. 




Figure 2. MIPS 24 urn (a) and 70 urn (b) images (the latter 
smoothed using a box size of 2 pixels) of PKS 0043-42. Panels are 
110" side. 
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Figure 3. SED of PKS 0043-42 including the NIR Ks-band upper 
limit, the IRS observed spectrum scaled to the MIPS 24 um point, 
the MIPS 70 um upper limit, and the 15 and 2.3 GHz radio core 
fluxes. The latter indicate that non-thermal emission from the 
radio core does not contribute significantly to the MIR emission. 



ing IRS data for all the WLRGs in the 2Jy sample with red- 
shift 0.05<z<0.15 (excluding PKS 0034-01, which is affected 
by saturation problems), and is representati ve of WLRGs 
at sim ilar redshifts to PKS 0043-42 (see also iLeipski et alJ 
2009). At short wavelengths, this average WLRG spectrum 
shows a blue spectral slope, which corresponds to the tail of 
the photospheric emission from the stars in the host galaxy. 
At longer wavelengths, the continuum appears relatively 
flat, although rising slowly with wavelength. This shape is 
likely produced by a combination of non-the rmal emission 
and e xtended dust emission heated by stars |Leipski et al.l 
2009). As mentioned in Section[TJ there is evidence for strong 
synchrotron contamination at MIR and far-infrared wave- 
lengths in many WLRGs, based on their IR -to-radio SEDs 
(|Dicken et alj2008tlvan der Wolk et al.ll201fj| ). Finally, emis- 
sion lines such as [Ne II]A12.8, [Ne ffl]A15.5, and [O IV]A25.9 
urn are also detected in the median WLRG spectrum. 



also rule out significant contamination (by extended emis- 
sion from the radio lobes) of the IRS data (maximum slit 
width of 4.46"). 

The MIR spectrum of PKS 0043-42 (Figure [JJ shows 
a clear 9.7 [im silicate absorption feature and a spectral 
turnover at A > 20 u.m. These are the most convincing de- 
tections of torus signatures to date in a WLRG. We can 
quantify the 9.7 (im silicate feature in terms of its strength, 



defined as 



= ln(F ccmt )-m(F core ). We calculate F c 



using a cubic spline interpolation in logarithmic space over 
two intervals at shorter and longer wavelength than the re- 
gion covered by the silicate feature to fit the continuum 
and F eore allowing the peak wavelength to vary, following 
ISirockv et all (|200&t ). By doing this, we find t^ p =Q.2. The 
lack of strong emission lines in the IRS spectrum is also 
remarkable. 

All of these features, as well as the spectral shape, are 
markedly different from those of the median spectrum shown 
in the bottom panel of Figure [TJ which was constructed us- 



2.2 Chandra observations 

PKS 0043-42 was observed with Chandra on 2009 Feb 05 
(obsid 10319) as part of a systematic programme of Chandra 
observations of the complete steep-spectrum 2Jy sample in 
the redshift range 0.05 < z < 0.2 (details of all the observa- 
tions, including further results from the observation of PKS 
0043—42, will be presented by Mingo et al. in prep). The 
time on source was 18.4 ks. X-ray emission from the nucleus 
was clearly detected, at the level of ~ 150 counts in the 0.5- 
5.0 keV energy band, and we extracted a n X-ray spectrum 
of the nucleus in the manner described by lHardcastle et al.l 
(2006). This spectrum is shown in Figure [4] 

The spectru m is clearly different fro m what is usually 
seen in WLRGs l|Hardcastle et all l2006h . namely a single 
power law with little or no intrinsic absorption: a peak at 
high energies indicates the need for an absorbed component 
in the model. Therefore, we fitted the standard model for 
narrow-line radio galaxies, consisting of two power laws, one 
with absorption fixed at the Galactic value (Nh = 2.7 x 10 20 
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Figure 4. Chandra ACIS-S spectrum of PKS 0043-42. The model 
consists of local absorption, intrinsic absorption, two power-laws 
and a Gaussian to account for the Fe emission line. The presence 
of this emission line is associated with accreting structures, which 
are typically seen in NLRGs, but have not been unequivocally 
observed in WLRGs before. 



cm" ), and one with an additional free variable absorption. 
Because of the small number of counts in the spectrum, we 
fixed the photon indices of the u nabsorbed and a b sorbe d 
power laws to the values used by lHardcastle et all (2006), 
i.e. 2.0 and 1.7 respectively. A fit of this model showed strong 
residuals at 6-7 keV, which could not be removed by allow- 
ing the photon index of the absorbed power law to vary; a 
good fit could only be obtained with the addition of a nar- 
row Gaussian feature with E = 6.47to' 4 keV, consistent 
with being an Fe Ka line (indicative of the presence of an 
accretion disk). The best-fitting absorption column density 
in this model was OOlg'gg x cm_2 an d the unabsorbed 
2-10 keV luminosity of the heavily-absorbed power-law com- 
ponent was 1.831q '^2 x 10 43 er 8 s_1 - The high obscuring col- 
umn density and the clearly detected iron feature are both 
very commonly seen in NLRG s, but have not so far been 
seen m bona fide WLRGs (e.g. lHardcastle et alj[200sh . 

In addition, the X-ray luminosity lies in the region 
populated by NLRGs on a plot of absorbed X-ray lumi- 
nosity against low- fre quency radio luminosity (Fig. 3 of 
lHardcastle et al" . 2009), as well as on the correlation between 
the 15 u.m luminosity (5.8xl0 43 erg s -1 , as measured from 
the IRS s pectrum) and the unab sorbed X-ray luminosity 
(Fig. 10 of lHardcastle et alj|2009h . Thus, the X-ray results 
strongly support the picture in which PKS 0043—42, while 
spectroscopically classified as a WLRG, has a torus and ac- 
cretion disc more typical of a NLRG. 



3 SED MODELLING 

Our aim is to investigate whether or not the IR SED of PKS 
0043-42 can be reproduced by torus models. The IR range 
(and particularly the MIR; 5-30 um) is key for setting con- 
straints on the torus models, since the reprocessed radiation 
from the dust in the torus is re-emitte d in this range. Recent 
ground-based MIR im aging (see e.g. iPackham et al.l 120051 ; 
iRadomski et al.l l2008f) and interferometric observation s of 
nearby AGN (e.g. lJaffe et~aTll2004l ; iTristram et alj|2007h re- 



veal that the torus size is likely restricted to a few parsecs - 
much smaller than the resolution of the IRS spectra. Thus, 
large aperture data such as those employed here may po- 
tentially be contaminated by i) photospheric emission from 
stars, ii) extended dust emission heated by stars and/or 
the AGN, and iii) non-thermal emission. First, the optical 
spectrum and the lack of polycyclic aromatic hydrocarbons 
(PAHs) in the MIR suggest that star formation is not impor- 
tant in PKS 0043-42. Second, extended dust emission from 
the NLR cannot be a substantial source of contamination 
either if we consider t he low emission line luminosity mea- 
sured for this galaxy (|Tadhunter et all Il998l ; iDicken et al.l 
2009). Finally, in Section [2~T1 we have shown that the non- 
thermal emission from the radio core and the lobes does not 
contribute significantly to the MIR emission of this object. 
Thus, the dominant contribution to our IR data must be 
nuclear dust heated by the AGN. 



The clumpy dusty torus models of lNenkova et al] (|2002h 
hold that the dust surrounding the central engine of an 
AGN is distributed in clumps, instead of homogeneously fill- 
ing the torus volume. These clumps are distributed with 
a radial extent Y = R /Rd, where R and Rd are the 
outer and inner radius of the toroidal distribution, respec- 
tively. The clumpy database now contains 5xl0 6 models, 
calculated for a fine grid of model parameters. To take 
into account the inherent degeneracy between these param- 
eters when fitting observables, here we use an updated ver- 
sio n of the Bayesian inference tool Ba y esClu mpy, developed 
bv lAsensio Ramos fe Ramos Almeida! ([2009) . BayesClumpy 
applies interpolation methods that are able to derive models 
for different combinations of the clumpy model parameters 
even if they are not present in the original database. For 
an example of the use o f clumpy model fitting to IR SEDs 
using BayesClumpy, see I Ramos Almeida et al. (2009). The 
new version of BayesClumpy allows to fit photometric and 
spectroscopic data simultaneously. Thus, her e we use an in- 
terpola ted version of the clumpy models of iNenkova et all 
(2008a b) including the corrections for the previously erro- 
neous AGN scaling factor (see erratum by INenkova et al.l 
l2010h . 



The result of the SED fitting are the poste- 
rior/probability distributions for the six free parameters 
that describe the models (see Table Q] for a description of 
the parameters) and the vertical shift (Figure 0. In ad- 
dition, we can translate these results into corresponding 
spectra (Figure [BJ. Thus, the maximum-a-posteriori (MAP) 
spectral energy distribution (SED) represent the "best fit" 
to the data (solid line), i.e. the combination of parameters 
that maximizes the overall probability/posterior distribu- 
tion. The dashed line is the SED obtained using the median 
value of the probability distribution of each parameter, and 
the shaded region in Figure [6] indicates the range of models 
compatible with all the possible combinations of parameters 
within their 68% confidence intervals around the medians. 
Figure [6] shows how well the clumpy models fit the observa- 
tional data. Both the 20 [im turnover and the silicate feature 
in weak absorption are reproduced by the MAP and the me- 
dian models. 
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Figure 5. Probability distributions of the free parameters that describe the clumpy models resulting from the fit of PKS 0043-42. 
The vertical shift has been marginalized. Solid and dashed vertical lines represent the MAP and the median values of each distribution 
respectively, and the dotted vertical lines indicate the 68% confidence level for each parameter around the median. 



Parameter 


Abbreviation Interval 


Fitting results 






Median Mode 



Width of the angular distribution of clouds 


(7 


[15°, 75°] 


60°± 6 ° o 

13 

12±1 


68° 


Radial extent of the torus (Ro/R^) 


Y 


[5, 30] 


12 


Number of clouds along the radial equatorial direction 


N 


[1, 15] 


5±1 


5 


Power-law index of the radial density profile 


1 


[0, 3] 


0.3±g;l 


0.3 


Inclination angle of the torus 


i 


[0°, 90°] 


74° ±8° 


75° 


Optical depth per single cloud 


TV 


[10, 200] 


18±S 


18 



Ay produced by the torus along the LOS Ay OS . . . 84±y 1 mag 79 mag 

Table 1. Clumpy model parameters. Columns 1 and 2 give the parameter description and abbreviation used in the text. Column 3 
indicates the input ranges considered for the fit (i.e., the lower and upper limits of the uniform priors). Finally, columns 4 and 5 list the 
resulting median and mode values of resulting probability distributions resulting from the fit (Figure [5j ■ 



4 RESULTS 

The more information provided by the IR SED, the better 
the probability distributions are constrained. From Figure 
\5\ it is clear that the IR data of PKS 0043-42 constrain the 
six parameters well (two last columns of Table [TJ). Thus, 
the IR emission of this radio galaxy can be reproduced by a 
broad clumpy torus (torus width of <t=60°) with an average 
number of clouds No=5 along the radial equatorial direc- 
tion, a radial extent Y=12, and close to edge-on inclination 
(i=74°). The optical extinction produced by the torus along 
the line-of-sight (LOS) is A v os = 1.086 A^o r v e ( - (l - 90)2 /ct2) 
= 84 mag. 

The clumpy model fits yield the intrinsic bolomet- 
ric luminosity of the AGN by means of the vertical shift 
applied to match the observational data points: ^ii R = 
1.6±o;i x 10 44 erg s~ 1 . This value can be directly com- 
pared with the unabsorbed 2-10 keV bolometric luminos- 
ity measured from Chandra hard X-ray o bservations by us - 
ing a bolometric correction factor of 20 (|Elvis et al.lll994T ): 



L^, = 3.7xl0 44 erg s" 1 (see Section Thus, we find 

L* ol ~ 2xL^ ; Ii? , which is consistent within the uncertain- 
ties. In fact, these lum inosity values are typical of SLRGs 
i|Hardcastle et al.ll2009l ). and both the MIR (2 4 and 70 urn) 
and [O IIIJA5007 A luminosities reported in iDicken et alj 
( 2009) for PKS 0043-42 are at the upper end of the WLRG 
range in the 2Jy sample. 

The outer size of the torus (R — YRd) scales with the 
AGN bolometric luminosity, so assuming a dust sublimation 
temperature of 1500 K, R a = 0.4 Y (L^^/IO 45 ) 5 pc. 
Considering the median value of y=12 determined from our 
fit, we find that the torus in PKS 0043-42 has an outer extent 
of R ~ 2 pc. This value is consistent with interferometric 
observations of nearby AGN, which indica te that the torus 
emission only extend s out to 1-2 pc (see e.g. lJaffe et al.ll2004l ; 
iTristram et alj|2007h . 

Using the A v os value reported in Table [1] we can de- 
rive the column density usin g the Galactic dust -to-gas ratio 
(N|° s = 1.9 x 10 21 Af, os ; lBohlin et al.lll97Sl ). This gives 
a value of N| os = 1.6 ±g;? xlO 23 cm" 2 , which is simi- 
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Figure 6. Fit of the observed IRS spectrum (thick solid line span- 
ning 5-30 (im) and the SOFI Ks and MIPS 70 urn data. Dashed 
and solid lines correspond to the median and MAP models re- 
spectively. Shaded region indicates the range of models compat- 
ible with all the possible combinations of the parameters within 
their 68% confidence intervals around the medians. 

lar to that measured from the Chandra X-ray observations 

(^X-rays = ^ ± 0.S x 1Q 23 ^-2. Section . This indi _ 

cates that the columns of material implied in the X-ray ab- 
sorption are comparable to those inferred from our IR data. 
Thus, dust-free gas absorption (i.e. gas located inside Kd) 
does not appear to be important here. This result contrasts 
with the findin g s for S eyfert galaxies and QSOs reported by 
iMaiolino et al.l ()200ll ). who claimed that the absorption of 
the nuclear region in these objects, as derived from optical 
and IR broad lines, is generally much lower than that ob- 
tained from the gaseous column density deduced from the 
X-rays, if a standard Galactic dust-to-gas ratio is assumed. 



5 DISCUSSION AND CONCLUSIONS 

Based on our analysis, we conclude that the MIR emission 
of PKS 0043-42 can be accounted for by a clumpy dusty 
torus with the characteristics shown in Table [T] The derived 
model parameters are compatible with those obtained for 
Ty pe-2 AGN of similar AGN b olometric luminosities (see 
e.g. iRamos Almeida et alJliobH ). The presence of a torus is 
also confirmed by hard X-ray Chandra observations, which 
show clear evidence for a heavily absorbed AGN nucleus 
and a strong iron K-alpha line, which are commonly seen in 
SLRGs but not in WLRGs (Section Indeed, its AGN 
bolometric luminosity (L^,; = 3.7xl0 44 erg s _1 ) indicates 
the presence of a moderately luminous hidden AGN. It has 
been proposed that the properties of WLRGs, whose spectra 
lack prominent emission lines, are the result of the Bondi ac- 
cretion of hot gas, which prevents the formation of the clas- 
sical AGN structures such as the torus or the BLR. However, 
the presence of a compact dusty torus and accretion disk sig- 
natures in PKS 0043-42 provides evidence that this WLRG 
is fuelled by cold, rather than hot gas accretion. 

One question that arises from our result is: if there is a 
moderately luminous AGN heating a nearly edge-on dusty 
torus, why are we not detecting the strong, high equivalent 



width emission lines characteristic of NLRGs at either opti- 
cal or MIR wavelengths? It is important to emphasise that, 
although PKS 0043-42 falls close to the correlation between 
MIR and [O III] luminosity for 2Jy radio galaxies, its [O III] 
luminosity is a factor of ten lower than that of SLRGs a t 
similar redshifts and radio powers |Tadhunter et all ^998). 
leading to the low equivalent width of the emission lines at 
optical wavelengths. On the other hand, in the case of the 
MIR spectrum, the non-linearity o f the correlation between 
emission line and MIR luminosity l|Dicken et alj|2009l . I20T0I ) 
will naturally lead to the objects at the lower luminosity 
end of the correlation having smaller emission line equiva- 
lent widths than their higher luminosity counterparts. In the 
case of PKS 0043-42, this effect is exacerbated by the fact 
that it lies slightly above (~0.3 dex, or a factor of two) the 
L24p.m vs. Lyoin] correlation, potentially making the MIR 
fine structure lines even harder to detect against the thermal 
continuum. 

Note that the position of PKS 0043-42 above the Z/24/jm 
vs. Lyoiu] correlation can be explained if the dusty torus 
structure emitting the MIR continuum has a relatively large 
covering factor com pared with the NLR (see discussion in 
iDicken et al J 120091 ) . Alternatively, if the illuminating AGN 
is highly variable and has recently "switched on", the MIR 
continuum will be artificially enhanced relative to the emis- 
sion from the NLR, because the compact nuclear dust struc- 
ture (r~2 pc according to our modelling) will respond more 
quickly to the change in AGN illumination than the larger- 
scale NLR (r~0.1-l kpc). 

It is also interesting to consider the relationship be- 
tween PKS 0043-42 and the other WLRGs in the 2Jy sample. 
Along with the marked differences in the shape of its MIR 
continuum spectrum demonstrated in Figure[T]and discussed 
in Section 1231 PKS 0043-42 is 2 - 10 times m ore luminous 
at 24 urn than the other WLRGs at z < 0.2 (|Dicken et al.1 
2009). Moreover, for the majority of these low redshift WL- 
RGs, the overall MIR continuum shape and extrapolation of 
the radio core contiuum to MIR wavelengths suggest that a 
substantial fracti on of their MIR co ntinuum represents ei- 
ther synchrotron (|Dicken et al.|[2008h or ther mal dust emis- 
sion associated with starburst components l|Leipski et al.l 
2009). Therefore, the position of most of the WLRGs close to 
the Z/24fim vs. Liqiii] corre lation derived for th e 2Jy sample 
as a whole (see Figure 6 in lDicken et al.ll2009T ) may be mis- 
leading, since they would lie below the correlation if the syn- 
chrotron and the starburst contributions were subtracted. 
Indeed, there is no evidence for warm dust emission from a 
torus in the spectra of any of the z < 0.2 WLRGs in the 2Jy 
sample apart from PKS 0043-42. If present, the torus com- 
ponent is likely to be sustantially weaker in most WLRGs 
than it is in PKS 0043-42. 

Overall, while PKS 0043-42 is likely to represent a low 
luminosity counterpart of the SLRGs in terms of its nuclear 
structures and fuelling via cold accretion, such continuity in 
physical processes is unlikely to hold for many low-z WLRGs 
which show no evidence for warm dust emission associated 
with the torus. Therefore our results provide strong evidence 
that not all WLRGs are fuelled via the same mechanism. 
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